arXiv: 1508.06361vl [cond-mat.mtrl-sci] 26 Aug 2015 


Non-circular skyrmion and its anisotropic response in thin films of chiral magnets under tilted 

magnetic field 

Shi-Zeng Lin* and Avadh Saxena 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

(Dated: August 27, 2015) 

We study the equilibrium and dynamical properties of skyrmions in thin films of chiral magnets with oblique 
magnetic field. The shape of an individual skyrmion is non-circular and the skyrmion density decreases with the 
tilt angle from the normal of films. As a result, the interaction between two skyrmions depends on the relative 
angle between them in addition to their separation. The triangular lattice of skyrmions under a perpendicular 
magnetic field is distorted into a centered rectangular lattice for a tilted magnetic field. For a low skyrmion 
density, skyrmions form a chain like structure. The dynamical response of the non-circular skyrmions depends 
on the direction of external currents. 
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A magnetic skyrmion is a topologically protected spin tex¬ 
ture, which has been observed in magnets without inversion 
symmetry recently such as MnSi and FeGe.^“^ A skyrmion is 
characterized by a topological charge Q = ^ j dr^S • (dxS x 
= +1 with S(r) being a unit vector describing the direc¬ 
tion of the spin. The typical size of skyrmion is about 5 nm to 
100 nm and skyrmions form a triangular lattice. In bulk crys¬ 
tal, the skyrmion lattice is stabilized in a small region close 
to the critical temperature in the temperature - magnetic field 
phase diagram.^ Skyrmions are found to be much more stable 
in thin films.Skyrmions respond to various external stimuli, 
such as magnetic field, electric current and temperature gra¬ 
dient. One extremely attractive feature of skyrmions is that 
they can be depinned by a low current density of the order 
of 10^ A/m^, which is 5 to 6 orders of magnitude smaller 
than that for magnetic domain walls.Moreover, the con¬ 
duction electrons in a metal interact with skyrmions and ac¬ 
quire a Berry phase, which produces an emergent electromag¬ 
netic field acting on these electrons. This gives rise to the 
topological Hall effect which has been observed experimen¬ 
tally.^"^^ Skyrmions in insulators can be driven by a tempera¬ 
ture gradient^^"^"^ or electric field^^’^^. For their unique physi¬ 
cal properties, skyrmions are believed to be a prime candidate 
for the next generation spintronic devices. 

It is crucial to tailor the skyrmion structure to optimize the 
desired functionalities. For instance, to achieve high density 
memory utilizing skyrmions, one needs to have skyrmions 
with the size in the nanometer range. The size of the skyrmion 
can be controlled by spin anisotropy or external magnetic 
fields, while the density of skyrmions can be tuned by exter¬ 
nal magnetic fields. Even the chirality of the skyrmion can be 
tuned by the sign of the Dzyaloshinskii-Moriya (DM) inter¬ 
action through chemical substitution.^^ The skyrmion in these 
cases has circular shape and the response is isotropic. It is of 
fundamental interest and of relevance for applications whether 
there exist non-circular skyrmions with an anisotropic interac¬ 
tion between them. We note that such non-circular skyrmions 
were observed experimentally in strained crystals, where the 
DM vector becomes anisotropic.^^ Here we propose a simple 
way to stabilize non-circular skyrmions in thin films of chiral 
magnet by tilted magnetic fields. 

Skyrmions are more stable in thin films by suppressing the 


competing conical phase when a magnetic field is applied per¬ 
pendicular to the thin film. When the magnetic field is tilted 
from the normal, the skyrmion phase is suppressed, which 
causes the decrease of the topological Hall resistivity as was 
measured recently in experiments.^^ When the field is paral¬ 
lel to the film, the skyrmion phase is suppressed completely 
and the conical phase is stabilized. Meanwhile the skyrmion 
shape is distorted for a tilted magnetic field because the re¬ 
gion with spin parallel to the in-plane component of the field 
grows while the region with the opposite spin shrinks. Be¬ 
cause of the distortion of skyrmion shape, the pairwise in¬ 
teraction between two skyrmions also becomes anisotropic, 
i.e. the interaction energy depends on the relative angle be¬ 
tween the two skyrmions. The resulting skyrmion lattice is 
no longer a triangular lattice. For a low density of skyrmions, 
chains of skyrmions are stabilized because of the anisotropic 
interaction. At last, we will show that the dynamical re¬ 
sponse of skyrmions to an external current drive also becomes 
anisotropic. 

We consider the following Hamiltonian for a classical spin 
S with IS I = 1 in two dimensions (x-y plane)^^ 

'H = -jY,SrSj-D Yj (S, X • e^) - H • ^S„ (1) 

(ij) i,fi=x,y i 

where J is the exchange interaction between the nearest neigh¬ 
bor spins, D is the DM vector along the bond due to the break¬ 
ing of inversion symmetry, H is the external magnetic field 
with a tilt angle 0 from the z axis, i.e. = F/^sin0 and 
Hz = Hcos(p. Here with yu = x, y is a unit vector in the 
A and y direction, respectively. The magnetic dipolar interac¬ 
tion was neglected because its strength is much weaker than 
the interactions in Eq. (1). Equation (1) can reproduce satis¬ 
factorily the measured phase diagram in experiments. To ob¬ 
tain the skyrmion configuration at zero temperature T = 0, we 
anneal the system using the Landau-Lifshitz-Gilbert equation 
with the Slonczewski spin-transfer torque term^^ 

5,S,- = -rS,x(Heff + H)+aS,x5,S,-^7ex,S,x(SyxS,). (2) 

The effective field is Hgff = -dlHIdSi and H is the Gaussian 
noisy field. Here a is the Gilbert damping coefficient and y is 
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Figure 1. (Color online) Skyrmion becomes non-circular when the magnetic field in the x-z plane is tilted away from the normal of the film. 
The tilt angle 0 is: (a) 0 = 0°, (b) 0 = 40°, (c) 0 = 60° (d) 0 = 66°. The color denotes the spin component along the z direction, while 
the arrows represent the in-plane spin component. Inset in (a) is a schematic view of the setup. Here N = 0.8T)^// and the plotted area is 
20 X 20J^a^/D^, with a being the lattice constant of the spin system. 


the gyromagnetic ratio. To discuss the skyrmion dynamics, we 
have also introduced the spin current /ext to describe the adi¬ 
abatic spin-transfer torque between the conduction electrons 
and localized spins. 

Let us first consider a single skyrmion in the tilted magnetic 
field. The spin configuration of a skyrmion as a function of tilt 
angle 0 is presented in Fig. 1. For a perpendicular magnetic 
field, the skyrmion is centrosymmetric. As the magnetic field 
deviates from the normal of the film, the skyrmion elongates 
along the direction perpendicular to the in-plane component 
of the magnetic field due to the Zeeman interaction term. For 
a tilt angle 0 > 66°, a clear coexistence phase of the skyrmion 
and the conical state can be seen. Such a skyrmion remains 
metastable even for 0 = 90° when we tilt the field continu¬ 
ously towards the in-plane direction, because of the topologi¬ 
cal protection. 

For a non-circular skyrmion, the interaction between 
skyrmions becomes anisotropic. To calculate the pairwise in- 
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Figure 2. (Color online) Pair interaction U(Ri - R 2 ) between two 
skyrmions at different tilt angle of the magnetic field: (a) 0 = 0°, 
(b) 0 = 45° and 0 = 60°. Here li and I 2 are two vectors along the 
principal axes connecting the origin and the energy contour at C/ = 
-0.7974/)^//. As indicated by the ratio/ 2//1 > 1, the pair interaction 
becomes anisotropic when the magnetic field is tilted. Here H = 
O.^D^IJ, D - 0.2J, / = 0 and = Ja/D. We restrict the pair 
interaction in the region |Ri - R 2 I > 3/a//) because the method to 
pin skyrmions by fixing spins in the skyrmion core becomes invalid 
when they are close to each other. 


teraction, let us first define the center of mass of a skyrmion 
as its topological center 

R = / J dr^S • {d,S X dyS)r. (3) 

Here we have used the continuum approximation, S/ ^ S(r), 
which is valid when D/J 1. The interaction energy be¬ 
tween two skyrmions, U(Ri - R 2 ), depends on Ri - R 2 , in 
contrast to the |Ri - R 2 I dependence for circular skyrmions. 
We then calculate U (Ri - R 2 ) by fixing spins of a skyrmion 
around its topological center |r - Ri | < in order to pin the 
skyrmion at a desired position. The results for U(Ri - R2) 
are shown in Fig. 2. The interaction is mediated by ex¬ 
change of magnons between two skyrmions and is repul¬ 
sive. For circular skyrmions at a perpendicular magnetic 
field (0 = 0) the interaction is isotropic and is given by 
U(Ri - R 2 ) = U(\Ri - R 2 I) ~ i^o(|Ri - R 2 |/^) for a large sep¬ 
aration |Ri - R 2 I » JajD. Here the length scale ^ is related 
to the magnon gap, a is the lattice constant of the spin system 
and Kq{x) is the modified Bessel function. For a tilted mag¬ 
netic field, the interaction becomes anisotropic. The repulsion 
between skyrmions in the direction where they are elongated 
is stronger than that in the other directions. This can be best 
seen by looking at h/k, where k (k) is the vector along the 
principal axes connecting the origin and the energy contour at 
U = -0.7974/)^//. Here magnetic field has the component in 
the k direction. The ratio k/k increases with the tilt angle 0, 
meaning that the interaction becomes more anisotropic for a 
larger tilt angle. 

Such an anisotropic interaction between skyrmions has 
profound consequences for the equilibrium configuration of 
skyrmion lattice. The triangular lattice is distorted into a cen¬ 
tered rectangular lattice (space group p6mm to c2mm transi¬ 
tion), i.e. the skyrmion lattice constant in the direction with 
a stronger repulsion becomes larger than that in the other 
directions. This can be seen from the spin structure factor 
{Sz(-(i)S displayed in Fig. 3 (a). Still we have six domi¬ 
nant Bragg peaks and higher order peaks for the skyrmion lat¬ 
tice. The ratio k/k is larger than 1 because of the distortion 
due to the anisotropic interaction. In the coexisting phase of 
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Figure 3. (Color online) (a) Spin structure factor (S'^(-q)S' 2 (q)> plotted in the logarithmic scale at // = 0.336/, D - 0.7265/, T - 0.05/ 
and 0 = 45° obtained by Monte Carlo simulations of model Eq. (1) with L = 60a. The skyrmions form a centered rectangular lattice with 
l\ - 0.0642 X Inja and I 2 - 0.09016 x Inja. (b) A skyrmion chain state obtained by numerical annealing of model Eq. (2) to / = 0 at 
H - 0.1 I J, D - 0.2/, 0 = 48° and L = 60 JalD. The color denotes the spin component along z direction, while the arrows represent the 
in-plane spin component. 


skyrmions and ferromagnetic state with a low skyrmion den¬ 
sity, as shown in Fig. 3 (b), skyrmions arrange themselves 
into a chain along the direction with a weaker repulsion. We 
note that such skyrmion chains have also been observed in 
nano wires due to the geometry confinement.^"^ 



Figure 4. (Color online) (a-c) Temperature-magnetic field (T-H) 
phase diagram of the model Eq. (1) at different tilt angles of the 
applied magnetic field obtained by Monte Carlo simulations. Here 
D = 0.1265J and the system size is L = 60a. The phase boundary is 
obtained by checking the spin structure factor and spin configuration, 
(d) Skyrmion density as a function of tilt angle obtained by numeri¬ 
cal annealing of the model from the paramagnetic state to / = 0 at 
//= 0.7DV/and D = 0.2/. 


Figures 4 (a-c) summarize the equilibrium phase diagram 
of model Eq. (1) at different tilt angles 0. For a low tilt an¬ 
gle, we have a magnetic spiral with the ordering wavevec- 
tor q along the [11] direction at low magnetic fields. For 
intermediate magnetic fields, the skyrmion lattice was stabi¬ 
lized. The skyrmion lattice is distorted into a centered rect¬ 
angular lattice for an oblique magnetic field. At high fields, 
we have a field-polarized ferromagnetic state. The skyrmion 
lattice phase shrinks with the tilt angle, and it disappears at 
0 60°. In this case the spiral with q along the [11] di¬ 

rection transits into the conical phase with q parallel to the 
in-plane component of the magnetic field direction. (Here it 
is along the v direction.) The skyrmion density as a function 
of the tilt angle 0 at T = 0 obtained by numerical annealing 
of Eq. (2) is presented in Fig. 4(d). The skyrmion density de¬ 
creases and vanishes completely at 0 60°. The threshold tilt 

angle where skyrmion density vanishes decreases with tem¬ 
perature. The decrease of the skyrmion density as evidenced 
from the topological Hall resistivity as a function of tilt angle 
was measured in Mno. 96 Feo. 04 Si thin films recently in Ref. 21, 
where the topological Hall resistivity vanishes at 0 40° at 

r = 20 K. 

Finally, let us discuss the equation of motion for the non- 
centrosymmetric skyrmion. It is more convenient to adopt the 
continuum approximation here. We follow Thiele’s collective 
coordinate approach^^ by treating a skyrmion as a rigid object, 
i.e. S(r, t) = S(r - \t). In this rigid skyrmion approximation, 
dtS(r, t) ^ -(v • V)S and S x Hgff = 0. After properly inte¬ 
grating out the internal degrees of freedom for skyrmions, we 
obtain the equation of motion for skyrmions as particles^^"^^ 
from Eq. (2) 


^ij (‘^ext,; + — 0- 


(4) 
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Figure 5. (Color online) Skyrmion form factors defined in Eq. (6) 
and the ratio of the Hall angle defined in Eqs. (7) and (8) for external 
currents in the x andy directions. Here a - Q.\,H - O.SD^//, T - Q 
and D - 0.2/. 

Here /, j = x, y and summation over repeated indices is as¬ 
sumed. The tensor G and the form factor tensor t] are given 
by 

= ~o)’ 

riij = T J dr^diS ■ djS. (6) 

For a non-circular skyrmion, ^ rjyy and r]j(y ^ 0. We com¬ 
pute numerically rjij for skyrmions in Fig. 1 at different 0, 
and the results are shown in Fig. 5. While rj^y ~ 0 and 
is almost independent on 0, rjyy increases rapidly with 0 when 
the conical phase starts to appear. Note that y is the direction 
along which skyrmions are elongated. For a current in the v 


direction, the skyrmion acquires a velocity in the y direction 
because of the damping. The Hall angle Ox of the skyrmion 
motion is defined as 

tan(0;c) = Vylv„ - + arjyy). (7) 

For a current in the y direction, the Hall angle is given by 

tan(0j,) = -vjvy = -arjyyiil - (8) 

As shown in Fig. 5, the response of skyrmions to current be¬ 
comes anisotropic, i.e. the Hall angle depends on the direction 
of the current, because of the non-centrosymmetric nature of 
the skyrmion in tilted magnetic fields. 

To summarize, we have studied the equilibrium phase 
and dynamics of skyrmions with an oblique magnetic field. 
When the magnetic fields are tilted away from the nor¬ 
mal of magnetic films, the skyrmion phase is less favor¬ 
able and the skyrmion density decreases. Meanwhile, the 
shape of a skyrmion becomes non-centrosymmetric, which re¬ 
sults in an anisotropic interaction between skyrmions. This 
anisotropic interaction stabilizes a centered rectangular lattice 
of skyrmions. In the low skyrmion density regime, a chain of 
skyrmions can be stabilized in the ferromagnetic background. 
Note that skyrmion chains have been observed in a confined 
geometry recently^ . The Hall angle of skyrmion motion de¬ 
pends on the direction of the current relative to the magnetic 
field direction. The predicted non-circular skyrmion and the 
resulting skyrmion configuration can be checked by imaging 
methods, such as Lorentz transmission electron microscopy 
or magnetic force microscopy. 
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